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Abstract 

We present the production cross-section of 7-rays based on data of p-p colli- 
sions at the Large Hadron Collider (LHC), revising the previous semi-empirical 
formula mainly for 1) the inelastic cross-section in p-p collisions, erpp(i?o), and 
2) the inclusive 7-ray spectrum in the forward region, <jpp^ 7 (E , E 7 ). We find 
that the previous cross-section gives a significantly softer spectrum than found 
in the data of LHC. In this paper, we focus our interest mainly upon the LHC 
forward (LHCf) experiment, giving 7-ray spectra in the very forward region 
with the pseudo-rapidity rj* > 8.8 in the center of mass system (CMS), which 
have not been reported so far. We also give the pseudo-rapidity distribution of 
charged hadrons with -3 < rj* < 3 obtained by ALICE and TOTEM experi- 
ments, both with LHC. We find that the revised cross-section reproduces quite 
well the accelerator data over the wide energy range from GeV to 30PeV for 
projectile protons, corresponding approximately to 100 McV to 3PcV for sec- 
ondary 7-rays. The production cross-section of 7-rays produced in the forward 
region is essential for the study of 7-ray astronomy, while not important are 
those produced in the central region in CMS, and of much less importance in 
the backward. We discuss also the average transverse momentum of 7-rays, p^, 
and the average inelasticity transferred to 7-rays, k*, obtaining that the former 
increases very slowly with p t = 100 - 220MeV/c_for E = 1 GeV - 26PeV, 
and the latter is almost independent of E , with k* 1/6, while we can not 
exclude the possibility of a small increase of k*. 
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1. Introduction. 

The high energy 7-ray observations open a new window not only for the 
astronmy, but also for the particle physics, cosmology, which must bring us crit- 
ical information and hints for the understanding of current questions such as 
nature of the dark matter, mechanism of the 7-ray burst, origin of the highest 
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cosmic-ray (CR), and so forth. In recent years, remarkably developed are tech- 
niques in both on-board and ground-based telescopes, Fermi, H.E.S.S., MAGIC, 
VERITAS, which have reported exciting results for the sky- map of GeV-TeV 
7-rays in space, still working continuously, see Q for review. 

A new program called CTA (Cherenkov Telescope Array) is further progress- 
ing in the form of the international collaboration consisting of many scientists 
from Europe, USA and Japan 28]. It covers very wide energy range of 7-rays 
with several tens GeV to more than 10 TeV with much higher sensitivity than 
ever achieved, the full operation of which is scheduled around 2020. If the CTA 
program is operated as scheduled, we expect to detect TeV 7-ray sources more 
than thousands, comparable with those currently observed in GeV region, which 
must bring us surely new aspects for the understanding of the universe. 

Under these situations, it is quite desireable to have a reliable production 
cross-section of hadron-induced 7-rays with TeV energy or more, while those of 
electron-induced ones are, needless to say, well established on the basis of QED. 
Stecker 0] presented the production cross-section of 7r°, o-pp^ w o(£ , , E v o) with 
E„a < 100 GeV (Eq < 1 TeV) in the form of very useful parameterizations in 
1973, and later one of the authors (T. S.) revised it in the previous paper [l[ 
(hereafter Paper I) , applicable for the higher energy region Eq > 1 TeV. Nev- 
ertheless, the reliability is not always satisfactory, particularly in the forward 
region for E 1 > 1 TeV (Eq > 10 TeV) because of the limitation in the accelerator 
data, having no data on the energy spectra in the forward region, while Chacal- 
taya emulsion chamber (EC) group [20, 33] have given them using CR-beams, 
but poor in data quality. 

Fortunatelly, the LHC experiment started the operations in early 2010, 
among which the LHCf group measured 7-ray spectra in the very forward region 
with the pseudo-rapidity rf > 8.8 in the CMS [4[. The principal purpose of the 
LHCf is to find the best simulation code in hadron-interaction models, which 
plays a key role in the study of the extensive air shower (EAS) phenomena 
induced by the ultra-high energy CRs, affecting directly the estimation of the 
primary CR energy as well as its composition. While none of models currently 
available reproduce satisfactorily the LHCf data according to their preliminary 
studies HI], they will report a revised model soon based on further analyses as 
well as on coming data at ^/s = 14 TeV, the run of which is scheduled in 2014. 

Alternatively, the LHCf data also give us crucial information for the study of 
hadron-induced 7-rays in galactic environments, both interstellar medium (ISM) 
and the source of CRs, typically the supernova remnant (SNR). Particularly 
interesting is the production cross-section of 7-rays with 100 TeV or more in 
the laboratory system (LS), corresponding to the knee energy around PeV. We 
expect that the knee problem in close connection with the acceleration limit of 
CRs in the SNR may be solved by the CTA program through the observations 
of ultra-high energy 7-rays with J5 7 > 100 TeV, possibly much more clearly than 
those of hadronic components. Note that it is still not cleared albeit so many 
years have passed since the discovery of the knee [3(| . This is mainly due to the 
difficulty in observing the latter components around the knee energy by both 
direct (balloon and/or satellite) and the indirect (EAS) methods, which have 
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Table 1: Numerical values of [<xo, eo, k] appearing in Eqs. (1) and (2) for three cases of target 
mass, 1 < Ar£ < 4, 4 < A>p < 38, and 38 < Arp, where erg is in units of millibarn (mb), and 
eo in GeV. 



1 < Ap < 4 4 < A T < 38 38 < A T 
cq 20.0 2^3 18J3 
e : 0.93 0.92 0.92 
k: 2.23 2.40 2.40 



inevitable weaknesses in statistics for the former, and in the uncertainty of the 
composition for the latter. 

Now the accelerator data on c7pp_j. 7 (_Eo , E 7 ) being available over extremely 
wide energy ranges, Eo = 1 GeV ~ 30 PeV, it is an easy task to find empirically 
the cross-section by interpolating them without asking for complicated models 
in the meson production. While many phenomenological models with QCD 
have been applied for the current simulation codes [29(, one should keep in 
mind the fact that even the inelastic collision cross-section based on QCD, 
o-pp{E ), is not yet definitely established, and much less successful are those 
for the multiple meson production, app- > . 7T (EQ 1 E v ). In the present paper, we 
give CTpp_,. 7 (-Eo, Ey) with simple parameterization based on the experimental 
data with LHC, revising slightly the previous one, which will be quite useful for 
the study of future 7-ray astronomy even around 100 TeV or more. 



2. Cross-sections 

2.1. Inelastic collision cross-section 

In Paper I, we gave the empirical formula for the inelastic cross-section, 
(jpp(iSo), m P"P collision which covers the wide energy range from the threshold 
energy of pion production {yfs ~ 2 GeV) to the FNAL energy (y 7 ? = 1.8 TeV). 
Now we have the LHC data (ATLAS [5], ALICE[23], TOTEM[24]) on o- pp at 
^/i = 7 TeV. The LHC energy currently available corresponds to approximately 
26 PeV proton in the LS, hign enough even for 100 TeV 7-ray astronomy. 

Based on the LHC data, we assume a following empirical form for crpp, 



cTpp(£ ) = £o(£o)[i+ ys/^J, (1) 

with Z (E ) = o-o(c/vr(l-e- E °^°y (2) 

where Eq and v are the kinetic energy and velocity of the projectile proton in LS 
respectively, ^/so — 156 GeV, and see the first column (1 < Arr < 4) of Table 1 
for (Tq, eo, and k. Practically, Sq ps ctq for Eq > 10 GeV. 
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In Fig. 1, we demonstrate Opp against y/s including the LHC data, together 
with the previous empirical curve (dasshed curve), and the present one (solid 
curve) given by Eq. (1) with the numerical values summarized in the first column 
of Table 1. One sees that the previous one gives significantly over-estimation in 
the LHC energy region. 

For the inelastic collision between nucleus i and nucleus j, we give the cross- 
section based on the optical model [6, 7, 8], 

cry = 7r(a 4 2 + a 2 ) [In Xij + E x (xij ) + 7 E ] , (3) 

with 



v - A A gNN (A\ 
X«-A % A ^ + ,4) 



where Ai (Aj) is the mass number of the projectile (target) nucleus, 7g the 
Euler constant (=0.5772), E\{\) the exponential integral function, and a; (a,j) 
is related to the nuclear root- mean-square radius of the nucleus i (j), see [7| for 
the explicit value of a% (dj). Here cr^j^ is given by Eq. (1), but numerical values 
of the parameters appearing there depend slightly on the target nucleus A^ as 
presented in Table 1. 

2.2. Production cross-section of "{-rays 

Now in the following discussions, we use the natural units with c — 1 (the 
speed of light) unless otherwise mentioned specifically, and the asterisk attached 
to variables denotes those in the CMS, or else those in the LS. 

In Paper I, we assumed that the distributions of the energy and transverse 
momentum in the multiple meson production are both given by the exponential 
function, exp[— (E*/T n + p^/p )} (E*: 7-ray energy, nr. its transverse momen- 



tum), which are expected from the fire-ball picture 12] . On the other hand in 
the QCD picture, the energy distribution is given by the algebraical function, 
(I — x*) m , expected from the quark dimensional counting approach [l3j], where 
x* is the light cone momentum fraction, and m relates to the number of quarks 
actively involved in the collision with m = 3 ~ 6 practically. 

While both types of the distribution are equivalent in the central region 
around x* ~ 0, we find that the algebraical- type reproduces rather well the 
LHCf data in the exteremely forward region with x* w 1 (see Section 3.3). 
So in the present paper, we assume a following distribution function for the 
invariant production cross-section of 7-rays, modifying slightly the functional 



Figure 1: Inelastic cross-section <rpp in proton-proton collision as a function of the center of 
mass energy y/s ([9-11], [5], [23], [24]) with two empirical curves, the present one (solid curve) 
and the previous one (dasshed curve). 



5 



form used in Paper I, 



^ PP 7 d 3 p; 4vrT c 2 X* ' l ° j 

with 

x*=E*/T c ; Tc = ^/2 = Mp 7 c, (6) 

X* = x* + (Pt/Po; P t = E;sm6*, (7) 

where 7c is the Lorentz factor of CMS against the LS, Mp the proton mass, 
and 6*c is the normalization constant. 

In Eq. (5), we introduce four parameters, [m, C, jV 7 ,p ], where m corresponds 
to the softness of the energy-spectrum, £ to the correlation strength between 
E* and p^, and 7V 7 is of course the multiplicity of 7's, and po links to the 
average transverse momentum p^ (see Eq. [15]). In the present work, however, 
we fix the former two, with m = 4 and £ = 0.02, expected from the preparatory 
calculations [3lj , and focus our principal work on the determination of the latter 
two, [N 7 ,p ], (practically [iV 7 ,p^], see the beginning of Section 3) by the least 
square method in fitting with the experimental data. 

In order to make following discussions easy to understand, we introduce a 
parameter tc, corresponding to to defined by Paper I, 

t c = T c /po = v^/>o, (8) 

and 9c is given byl f 1 F {t*) „ . 

Fe(r)= x\l°- x) 4 e- TX dx, (£ = 0,1). (10) 
with Jo 

In Table 2 we give the explicit forms of Fi(t) (£ = 0,1), together with those 
of the approximation around r w 0, and demonstrate the numerical value of 0c 
against tc for several choices of £ in Fig. 2. 



Table 2: Explicit forms of -Fo(t) and Fi(t), and their approximate expressions for r 0. 
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Figure 2: Normalization constant 6>c against Tc for several choices of £ . 



2.5. Energy- angular distribution 

Eq. (5) is rewritten with m = 4 as 

47tT c 1 + Ctq ' 

O^ffy^re^erffifefefefefpmpa" 8011 with experimental data in the 
Remembering the invariant phase space, E*dE*dQ* = E 1 dE 1 dQ, the above 
equation gives the distribution function in LS as 

d'N, _N 1 e c pl{i-xr e Y c _ Tex ^ (12) 



dE^dtt 2ttM p r e + (t 
with x — E 7 /Eo, and 

T 9 = 2( 7 2-l)(Mp/p o )sin0, (13) 

r0=2( 7 2-l)(l-/3 c cos0), (14) 
note rg w 1 + 7^0 2 for E ~> Mp in the forward region, leading to i w x*. 

2.4- Average transverse momentum and inelasticity 

In this subsection, we give the average transverse momentum of 7-rays using 
Eq. (11) (or Eq. [5]), which is a critical parameter in the multiple meson pro- 
duction, with a quite stable value, say 150MeV/c, almost independent of the 
interaction energy. It is immediately given by 

^ = 9 c f\* e ^\d{co S e*). (15) 
Po Jo 1 + Qt 6 

The total energy flow transferred to 7-rays in the CMS, EE*, is similarly 
obtained from Eq. (11), 

ee; e; Fl (rS) 



^Vc J c Jo l + C T e 

where E* is the average energy of 7-rays in the CMS. In Fig. 3, we present p^/Po 
and E*/T c simultaneously against tq for several choices of (, corresponding to 
Fig. 2. 

Defining the average inelasticity transferred to 7 -rays in the CMS, 

EE* 

*• " T^W (17) 
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we have 

*; = W^« (i 8) 

2 Jo 1 + C T e 

where we use the approximation y/s >• 2Mp for the practical purpose. Assuming 
k^ Q ^ = 1/2 (total inelasticity) and p n a = 1/3 (isospin symmetry), we expect 
k* — 1/6, which is discussed again in Section 4. 

2.5. Pseudo -rapidity distribution 

Practically in the accelerator data, we often use a variable of pseudo-rapidity 
n* instead of 9*, defined by if — — lntan#*/2. The pseudo-rapidity distribution 
is immediately given by, after integrating with respect to i* in Eq. (11), 

dN 1 _ N,6 C Fq(t*) _ _ r c 
dr]* 2 cosh 2 n* 1 + {t* ' ^ cosh?]*' 

It is also important to see the energy flow, d(EE*), transferred to 7-rays 
within (n* , n* + dn* ) , which is given by 

d(s E ;) _ N,e c F^r;) 

T c dn* 2 cosh 2 n* 1 + (t* ' 
see T T ft bie i( 2 r |^r fl ^T^^^i( m ). L s within (77*, rf + dn* ) is easily given by 



djEEj) _ coshfa* + n c )d{ZE*) 
dn* cosh n* dn* 



(21) 



fmiis 1 ^elsSrlSble+rel^A b efwteh^hi^fyefllr^fW^i.^^cf'ftnete^nQHi 
CMS, SE 1 = ^ c SE* r after integrating both sides of Eq. (21) over n* (sec Eq. 
[16] for the explicit form of EE*) because of the forward-backward symmetry 
in the CMS, leading to fc 7 w k*. 

In Fig. 4, we demonstrate three kinds of density distribution against n* 
simultaneously, p*^ (multiplicity), pj, (energy flow in CMS), and (energy 
flow in LS), for E = 10 11 , 10 14 , 10 17 eV, where the vertical axis is normalized 
to unity after integrating over n* , 

_ dN, ( d(EE*) d(EE, ( ) 



Figure 3: Numerical values of Pt/PO an d E*/Tc against tc for several choices of f, where the 
former corresponds to the left axis, and the latter to the right axis respectively. 
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Figure 4: Three kinds of density distribution (see Eq. [22]) against the pseudo-rapidity rf 
in the differential form, a) multiplicity, b) energy flow in CMS, and c) energy flow in LS, for 
three energies, E = 10 11 , 10 14 , and 10 17 eV. 



Figure 5: Three kinds of density distribution against the pseudo-rapidity rj* in the integral 
form, corresponding to those in the differential form presented in Fig. 4. 



On the other hand, we present the integral forms for three kinds of density, 
-P^f(> ? 7*), -Pg(>?7*), an d -Pg(>?7*) in Fig. 5, corresponding to Fig. 4. From these 
figures, one finds that the energy flow in LS comes from mostly those produced 
in the very forward region in CMS, particularly for higher energy. This means 
that essential is the production energy spectrum of 7-rays in the very forward 
region in CMS, contrarily not important in the central region, and of much less 
importance in the backward. 

While the present paper is focussed upon the 7-ray component produced 
by 7r°-decay, the pseudo-rapidity distributions of charged hadrons, dN^/dr/*, 
have been extensively studied with accelerator experiments, closely related to 
dN 7 /dri* given by Eq. (19). So we present the kinematical relation between 
them rather in detail in Appendix A, which is given by, assuming 7V 7 w ^ cn i 



dN. 



ch 



drj* 



l + ^V#^tanhV 
2 \ Po J F (t,.) 



dN, 



7 



drj* ' 



(23) 



where m T o is the mass of 7r°, and see Eq. (A12) in the appendix for A(r)*). 

3. Comparison with experimental data 

In this section, we determine two parameters appearing in Eq. (5), [Nj,po], 
while the other two, [to, £], arc fixed to [4, 0.02] as mentioned in Section 2.2. 
Now the average transverse momentum p^ is quite stable, and well established 
in both accelerator and CR experiments with 150 <~ 200 MeV/c. So we use the 
parameter p^ in place of po, which is given by Eq. (15), 

f-f^^^/f^^ (24 > 

Po Jo 1 + C T e / Jo 1 + Q T e 

right-hand side of Eq. (24). In practice, p is easily obtained by the iteration 
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Figure 6: Energy spectrum of 7-rays at Eq = 0.97 GeV ll jj . Empirical curves are given by 
Eq. (12) after integrating over cos0. 



Figure 7: Energy spectrum of 7-rays at Eq = 23.1 GeV [l5ll for different emission angles. 
Empirical curves are given by Eq. (12). 



method with the initial value of 200MeV/c for the set of [-BoiPtl ( or PciPtD' 
since Pq is also quite stable with 150-250 MeV/c for = 100-230 MeV/c (see 
Figs. 3 and 17). 

Explicit values of [-A^,^] are presented in each figure appearing in the fol- 
lowing subsections (see also Figs. 16 and 17, and Table 3), which are obtained 
by fiting the experimental data with the present empirical curve. 

3.1. The low energy region (Eq = 1 — 300 GeV) 

From Eq. (12), we can obtain easily the energy distribution in LS integrating 
over cos 9. In Fig. 6, we compare the empirical one thus obtained with the 
data at Eq = 0.97 GeV given by Bugg et al. [l4j3, and find that they are well 
reproduced with the numerical values of N 7 and presented in the figure. 

Fidecaro et al. [15| gave the production cross-section of 7-rays for different 
emission angles at Eq = 23.1 GeV in the LS, which is presented in Fig. 7 together 
with curves expected from Eq. (12). One finds the agreement is excellent for all 
emission angles. 

In Fig. 8, we show the energy distribution with use of the Feynman variable, 
Xp = x*cos9*, at two energies, Eq = 11.5, 203.7GeV [16, 17]. Empirical 
curves are obtained by replacing x* with xpsec9* in Eq. (11), where we must 
take care of the integral range for cos 6** with xp < cos 6* < 1. One finds the 
empirical ones reproduce well both data. In Fig. 9, we present the pseudo- 
rapidity distributions at Eq = 11.5, 203.7, and 299.1 GeV [16, 17, 18], together 
with the empirical ones obtained by Eq. (19). Our numerical curves are again 
in nice coincidence with the data. 

3.2. The high energy region (E = 0.5 - 200 TeV) 

In this subsection, we compare the experimental data in TeV region with 
our cross-section given by Eq. (5), which are obtained by ISR, FNAL, and 



1 Original data were given in the form of the production cross-section for n°, 
CT pp— Hr° C^O) E v o )• As they gave explicitly the number of events per 25MeV energy bin, 
we converted them into <Tpp_> 7 (Eo, E 7 ) by randomly sampling for the n° — > 27 decay in each 
energy bin. 
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Figure 8: Production cross-section of 7-rays with use of the Feynman scaling variable x 7 ^ at 
two energies, E = 11.5, 203.7 GeV [16, 17]. 



Figure 9: Pseudo-rapidity distributions at three energies, E$ = 11.5, 203.7, 299.1 GeV [16, 
17, 18]. Empirical curves are obtained by Eq. (19). 



Figure 10: Energy spectrum of 7-rays for several sets of the emission angle in the CMS at 
ISR energies [lst . Empirical curves are given by Eq. (11). 



the Chacaltaya EC with CR-beams. In Fig. 10, we present the ISR data |19f 
with y/s = 30.2, 44.7, and 52.7 GeV, each corresponding to E = 0.483, 1.06, 
1.48 TeV in the LS respectively, where empirical curves are obtained by Eq. (11). 
One might worry about some discrepancies appearing in the low energy region, 
E* < 1 GeV, but it is not so critical in the practice as mentioned in Section 2.5, 
namely important is only the high energy part in the forward region in CMS, 
see Fig. 5. 

As presented in Paper I, the Chacaltaya EC data [20j provide the fractional 
energy spectrum of 7-rays, / 7 = E 7 /EE y) in Eq = 30 — 200 TeV region. The 
relation between SE 7 and Eq is given by assuming the 7-ray inelasticity, fc 7 , 
while we have to take care of the bias-effect in EC experiments (k^ = 0.28), 
see Paper I for the detail. In Fig. 11, we show the / 7 -spectrum for both data 
and curves expected from Eq. (12) after integrating over the emission angle. 



Figure 11: Fractional energy spectrum of 7-rays obtained by Chacaltaya EC experiments 
[20j for three energy flow ranges, SE-y = 7-10, 10-20, and 20-50 TeV, each corresponding to 
the average energy for projectile proton, {Eg) = 29.2, 45.1, and 101.2 TeV respectively. 



Figure 12: Pseudo-rapidity distributions obtained by UA7 |22j |. and the Chacaltaya EC 
experiments |2(I . Empirical curves are given by Eq. (19). 
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Figure 13: Pseudo-rapidity distributions of charged hadrons obtained by UA5 [2ll . Empirical 
curves are given by Eq. (23). 



The agreement is quite well within the statistical error. 

We present the pseudo-rapidity distribution in Fig. 12 obtained by UA7 



22j, and EC data [20|, where we present the curve expected from Eq. (19) with 
y/s = 600 GeV. Again we find that the present curve reproduces nicely the data. 

In Fig. 13, we show the pseudo-rapidity distribution of charged hadrons 
obtained by UA5 [2l|, for three energies, ^=52.7, 200, and 546 GeV. The 
numerical curves are given by Eq. (23) with the assumption of = iV 7 , taking 
the 7T° — ► 27 decay into account. We find that they are in good agreement with 
the UA5 data. 



3.3. The LHC energy region (E = 0.4 - 30PeV) 

Now we compare our production cross-section with the LHC data most re- 
cently reported, while the final goal with y/s = 14 TeV will be available around 
2014. LHCf group [4] present recently the energy spectra of 7-rays in the very 
forward region, 77* > 8.8, at y/s = 7 TeV, corresponding to Eq = 26PeV in the 
LS. Let us apply our formula given by Eq. (11) for LHCf data, and estimate 

In Fig. 14, we give the energy spectra at y/s = 7 TeV, for two sets of \Arf , 
A(f>*}, (a) [77* > 10.94, 360°] and (b) [8.8K rf < 8.99, 20°], where two curves 
from our empirical cross-section are presented together. We find that they 
reproduce well the experimental data in spectral shape, but the absolute value 
of -ZV 7 (average photon yield) is of approximately 20% difference between them, 
53.4 for (a) and 68.0 for (b), while the latter is consistent with 71.5 expected 
from TOTEM with the pseudo-rapidity distribution of charged hadrons (see 
Fig. 15), assuming V 7 fs N c ^. 

In addition to the e nerg y spactra of 7-rays in the forward region obtained by 
the LHCf group, UA5 [HI, ALICE Q and TOTEM [24] present the pseudo- 
rapidity distribution of charged hadrons in the central region as shown in Fig. 
15, covering the energies y/s = 900 GeV, 2.36 TeV, and 7 TeV, where numeri- 
cal curves are obtained by Eq. (23) with the assumption of V 7 sa ^ch- *~*ne 



Figure 14: Energy spectrum of 7-rays obtained by LHCf at y/s = 7 TeV for two sets of [Ar)*, 
A<j>*], (a) [77* > 10.94, 360°] and (b) [8.8K 77* < 8.99, 20°], where separately presented are 
two detectors, Arml (open circle) and Arm2 (filled circle), each with the scintillation fiber 
and the silicon strip respectively. 
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Figure 15: Pseudo-rapidity distribution of charged hadrons obtained by UA5 [2lll . ALICE 
[13, and TOTEM [ii|_at y/s = 900 GeV, 2.36 TeV, and 7TeV, where curves are given by Eq. 
(23), assuming iV 7 -^ c h- 



finds that they are well in consistent with the experimental data, particularly 
interesting is that the concave shape around 77* ss is nicely reproduced. 



Table 3: Summary of [N-y, p^], where LHCf-1 corresponds to the data with [rj* > 10.94, 
Z\0*=36O°], and LHCf-2 to those with [8.8K rf < 8.99, A<j>* =20°], and "(ra)" (n = 
0,1,2,...) appearing in the column of Eo denotes "10"". Fig. 13 and Fig. 15 with aster- 
isk mark correspond to the data from the pseudo-rapidity distribution of charged hadrons, 
assuming iV 7 « N ch and p t (-y) fb PtM/ 2 - 
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Figure 16: The total production cross-section of 7-rays, jV 7 <Tpp, against the proton kinetic 
energy Eo, where those compiled by Stecker Q| (filled squares) are plotted together. The 
present curve (solid curve) is obtained by Eq. (1) for crpp and Eq. (25) for N-y respectively, 
where the previous one (dashed curve) is also presented. 



Figure 17: Average transverse momentum of 7-rays against Eo, where we present also those 
expected from charged pions (open circles) [37] with the assumption of ^(7) p^(ir)/2. 



3.4- The multiplicity and the average transverse momentum 

In Figs. 6-15, we present explicitly the numerical sets of [iV 7 ,p^] in the 
extremely wide energy ranges, Eq = 1 GeV ~ 26 PeV, which are summarized all 
together in Table 3. 

For the 7-ray astronomy, practically the most essential is the total production 
cross-section of 7-rays, iV 7 x opp, no matter how the emission-angle 9* (or the 
transverse momentum pi) appears in the functional form of the cross-section. 
After Stecker Q summarized it in 1973, we revised it in Paper I with the data 
covering TeV region but without LHC data. 

Let us present iV 7 <7pp against Eg in Fig. 16 with LHC data, using Eq. (1) 
for <7pp(i?o), where we give a solid curve obtained by the following empirical 
form for N~, 



N 7 (E ) =N E 



0.115 



exp(-0.47y£ ) 



(25a) 



N {E ) = 8.80 x 



l-exp(-0.15y.Eo) 



(25b) 



with Eq = Eq — 2m v in GeV, and presented together is a dashed curve from 
the previous parametrization [l( for Nj(Eo). We plot also iV 7 crpp (open circles) 
expected from the pseudo-rapidity distribution of charged hadrons, assuming 
7V 7 fa iV ch (see Figs. 13 and 15). 

One finds that the previous one gives significantly over-estimation in PeV 
region, and the present one reproduces nicely the experimental points in the 
very wide energy range, E — 1 GeV ~ 26 PeV. 
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Figure 18: The average inelasticity, k* transferred to 7-rays against Eq, where also plotted are 
those expected from the charged hadrons (filled squares) with the assumption of k* »i k*,/2. 



In Fig. 17, we give the average transverse momentum, p^, against Eq, to- 
gether with the empirical curve given by 



0.0286 



l-exp(-1.156V^^ 



(26) 



with 



po = m^c = 140MeV/c, 



while it is of little interest for the 7-ray astronomy, but important for the study 
of shower phenomena in the atmosphere. In Fig. 17 we plot also the half of 
the transverse momentum of the charged pions [37], assuming p^ (7) s=s p^(tt)/2. 
One finds that p^ increases slowly with E , as given by Eq. (26). 



4. Discussions 

In the present paper, interpolating experimental data nowadays covering the 
very wide energy range from GeV to 30 PeV, we have focussed our work on the 
construction of the semi-empirical formula for the inclusive production cross- 
section of 7-rays, crpp^. 7 (i?o, E-y), without asking for the cumbersome QCD- 
based models, and find that it reproduces excellently the machine data over the 
very wide energy ranges. 

The present simple parameterization in the formula should be compared to 
the simulation codes currently available in the CR community, which are usually 
very complicated, patching different models separately in low and high energy 
regions, and heavy in the sense that they are constructed so that all the com- 
ponents (n ± , 7T°, K ± , . . .) in both soft (small q 2 ) and hard (large g 2 ) processes 
simultaneously match with the accelerator data. So it is not an easy task for 
physicists other than a developer of simulation code to improve it freely, by 
contrast with the empirical formula, quite easy to touch the parameters appear- 
ing there. Of course one should keep in mind that our approach (interpolation 
method) is not valid for the study of extremely high energy shower phenomena 
in the atmosphere, say ^> 10 18 eV, where even the LHC can not cover, resulting 
in the need of some theoretical models in order to extrapolate the LHC data 
much higher, while the present parameterization is valid enough for the future 
7-ray astronomy, even up to PeV-7 observation. 

We have concentrated our interest upon two parameters, _/V 7 and p^, partic- 
ularly on the former. As mentioned often, the multiplicity N 7 plays an essential 
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role for the study of 7-ray astronomy, appearing always in the form of the total 
production cross-section, N 1 (E )xa-pp(E ) as presented in Fig. 16. 

Alternatively, 7V 7 is important also for the study of the shower phenomena 
in the atmosphere, appearing in the (total) inelasticity k^ ^ (sa ^£ot)' So let 
us present k* (see Eq. [18]) transferred to 7-rays against Eq in Fig. 18, where 
also plotted are those (filled squares) expected from charged hadrons (see Figs. 
13 and 15), assuming k* ss k\/2. An error-bar attached to each square comes 
from statistical ones to N 7 as presented in Fig. 16. One finds approximately k* 
« 1/6 as a whole, almost independent of Eq, while it is as small as 0.1-0.17 in 
the low energy region Eq < ITeV, and as large as 0.15-0.2 in the high energy 
region ^> 10 TeV, indicating a small increase as the energy gets higher. But we 
reserve the conclusion for future studies, either constant or the increase. 



It has been well-known that the attenuation length A for the intensity of CR 
hadronic components in the atmosphere is given hy X/ A = 1 — ((1 — fc^ot)^) 



25[, where A is the collision length, and (3 the index of the integral primary 



CR spectrum with ~ 1.8. Experimentally we have A/X 1=3 1.5, with for instance 
A w 70g/cm 2 and A ~ 100g/cm 2 [§7j|> leading to k\ ^ w 1/2 in the energy 
region Eq < 100 TeV, assuming the uniform distribution in fctat, wm le not yet 
clear in the air shower region Eq > 10 PeV. Anyway the present result is not 
inconsistent with the common understanding in the inelasticity expected from 
the attenuation of CRs in the atmosphere. 

Finally we address further two open problems in the present paper; (1) 
the nucleus effect of proton- nucleus (p-A) and/or nucleus- nucleus (A- A) colli- 
sions for the production cross-section of 7-rays in p-p collisions, and (2) the 
applicability of the present empirical cross-section for the galactic phenomena 
other than emissions of the hadron-induced 7-rays, particularly for those of the 
electron-positrons coming from tt— \i— c decays. 

First, for the problem (1), one should remember that the effective 7-rays 
produced by the nuclear interaction in the galactic environments (either in ISM 
or in SNR) are only those produced in the forward region in the CMS, while not 
important are those in the central and backward regions. The effect of the plural 
interactions inside the nucleus appears only in the latter regions. In fact, it has 
been experimentally well-known that the difference between those produced by 
p-p and those by p-A (A-A) collisions appears only in the latter regions, while 
they are well in coincidence with each other in the forward region. This is the 
reason why we present Fig. 5, stressing in the present paper how essential are 
the 7-rays produced in the forward region, in contrast not important in the 
central region, and of much less importance in the backward. 

Practically, of course, we need the production cross-section of 7-rays for p-A 
(A-A) collisions, as there exist additionally helium and heavier components in 
CRs (projectiles) as well as the helium gas in the ISM (targets), while unfortu- 
nately we have not yet a reliable model nowadays for the p-A (A-A) collisions. 
We have used the modified wounded- nucleon model of Gaisser and Schafer [34] in 
our past calculations [26, 27], where we mention that the uncertainty in the nu- 
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cleus effect is of the second order for the practice. In order to see the uncertainty, 
we have introduced so called the "enhancement factor" defined by e q — <7 a n/<Zpp, 
taking the energy dependence into account, where q a Y[ is the emissivity of 7-rays 
in the galaxy produced by all kinds of nuclear interactions with p-p, p-A and A- 
A, and qpp by those with p-p only. For instance, e q = 1.54 (Gaisser-Shafer [34]), 
1.50 (Cavallo-Gould [35]), 1.60 (Stephens-Badhwar [36]), and 1.53 (Shibata et al. 
[26]), indicating that the difference in the choice of nucleus interaction model is 
not so significant as compared to that in the choice of the propagation model. 
These results tell us also that the procedure in the calculation of the CR prop- 
agation becomes quite simple by the use of the enhancement factor e q . 

Second, for the problem (2), indeed we do not touch the intermediate me- 
son, 7T°, except the pseudo-rapidity distribution of charged hadrons (Figs. 13 
and 15), but the decay products 7's only, having no interest in the intermediate 
mesons, either via ir° or via heavier 77, rj', etc.. However, as long as focussing on 
the decay products such as electrons and/or neutrinos from muons, the present 
empirical form for 7's, erpp_j. 7 (-Eo, -E 7 ), is valid also for muons produced via 
pions, crpp_}. M (i?o, E^), while we have to take care of the mass difference be- 
tween photon and the muon. This is because both decays, ir° -^7 + 7 and 
7r ± — > /M^ + v (D), are isotropic two-body decays in the pion rest system, leading 
to the same kinematics in 7 and fi^ but different mass. Namely, we do not 
need the information of the intermediate pions also in the case of the muon 
production cross-section. This fact tells us that we have model-independently 
a kinematical relation between the emissivity of 7-rays and that of (secondary) 
electron-positrons in the galaxy, detail of which will be reported elsewhere in 
connection with the galactic electron-positron spectrum. 

In the near future, we will apply the present cross-section, erpp^ 7 (i?o, £7), 
for the observational data in TeV region currently available on both diffused 
7-rays and those from the source, while one of the authors (T. S.) have studied 



the former components in Fermi energy region, 100 MeV ~ 100 GeV [26j, [271 ] . 
using the old parameterization in Opp^ 7 (i?o, E-y)- 
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Appendix A: Kinematical relation between 7 and it angular-distributions 

The kinematical relation between 7 and 7r° energy-distributions was studied 
by Sternheimer [32l |. where he assumed that the opening angle, a*, of two 7's 
from the decay of high energy ir° is so small that the emission angle of 7 is 
approximately equal to that of ir° , 6* « 9* . But this approximation seems to 
be too rough to transform practically from the pseudo-rapidity distribution of 
7's to that of 7r°'s, namely dN^/drj* = 2dN 7r o/dr]* in his approximation. In this 
appendix we present a more realistic relation between them. 
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Let us consider a tt° produced by p-p collision with (E* , #* , 0* o ), each 
denoting the energy, emission angle, and the azimuthal angle in the CMS re- 
spectively, and disintegrate into two 7's with (E^,0l,(f>*) and (E^, , 0%, ^2)1 
where individual angles are those against the collision axis. We define further 
two sets of angles, ((9*, $1) and (6*2, $2); which are the emission angles and 
the zenith angles of two 7's against the moving direction of ir°. In this appendix 
we put c = 1 for the sake of simplicity. 

We have following relations in these variables 

e* 1+ e; = e; , (Aia) 

0* + OZ=a*, (Alb) 
E* sin 9* = E; sin 6> 2 * , (Ale) 

and 

cos 6{ = cos 9* cos 6»*o + sin 9* sin 6* cos . (A2) 
Integrating over for both sides of Eq. (A2), we have 

(cos 61) = cos 9* cos 6% , (A3) 

and hereafter we omit angle brackets, (•••), in the left-hand side for the sim- 
plicity, as we are not interested in the azimuthal component. 
Remembering a well-known relation in 7r°-decay [33J 

a ~ — , = — , A4 

^/EfE* y/EtiE'o-Ei) 

for E* ^> m^o, we obtain, from Eqs. (Ala)- (Ale), 



m^o E* 

& '~e*;\I^I lj 

namely 

cos 9{ ttl-u(E* ,El), 



with 



Now we have from Eq. (A3), putting 0* = 9* and E* = E*, 

coseio~{l+uj(E* 0l E;)}cos6;. (A5) 

Here we have to take care of the above expansion with respect to u>, which is 
based on the approximation with 9* 3> a*, while Sternheimer assumed a* ~ 
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32[ . This means the energy range of 7 is limitted within (E*_ , El ) for the fixed 
energy of ir°, E* , which are obtained from Eq. (A4) with 9* > a*, and given 

by 

E* ± = [1±B(E: o ,9* )}E: o /2, (A6) 

with 

b (e; (> ,81,) = v /i-(2m ffO /£;; o 0; o )2. (A7) 

Now, we define the energy-angular distribution function of tt° with i* = 

cos 6*o 



n n o (E*o,t*,o) 



dE* n dt* 



and thus the energy-angular distribution function of 7 is given by, putting t* = 
cos 0*, 

n~ / (E;,t*) = 2 1 f dtiosltio - \i + u(E; ,E;)}t;] n7r o(E* n0 ,tio), 

Je* P„o J-i > J 

(A8) 

and Tc is the maximum energy of 7r° given by Eq. (6) in the text. 
For E*„ 3> m^o (equivalently w < 1), we have with use of Eq. (A5) 



n„o (E* , t; ) « n„o (E* , t*) + u (E* ,E;)t;-^-n„o(E* ,t;), 



7 

and integrating over E* for both sides of Eq. (A8) in order to obtain the angular 
distribution, we obtain 

n^t*) « 2n w o(t;) + 2t; / ds; / —^. u (e^,e;) —n^(E; ,t;), (aq) 

Jo Je* P^o at 1 

where one has to take care of the kinematical constraints in (E* ,E*) given by 
Eqs. (A6) and (A7) for the practical integrations. 

Now we use the approximation given by Sternheimer for n^o appearing in 
the integrand in Eq. (A9), 



1 d 

n n o(E* ,t*) pa -~p n0 ^— nrr (E* ,t*), 



note that the second iteration for n^a is negligible as shown in Fig. 19, and next 
exchange the order of integrations for E* and E* a , taking care of the kinematical 
constraints mentioned before 

dE; / dE*o =^ / dE* / dE; 
Je* J e„, Je* 
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Figure 19: Correction rate for the Sternheimer approximation with a* ?s 0, corresponding 
to the second term in the square bracket in Eq. (All) for ^/s = 0.l, 0.2, 0.5, l.OTeV, where 
we assume po = 200MeV/c for two choices of £, 0.02 and 0.04. 



with 

E n * — 2r7j 7r o/#* = nijo/ tan~ 1 (e~' ? ). 



The second term in the right-hand side of Eq. (A9) is thus given by 

-^J^M^t^' 3 ^)^^'^ (A10) 

with 



x ' v ~ V \ x ) ' ,f ~ tan- 1 ^*)- 

Assuming 7V 7 = 2N w o , and substituting the explicit form of n 7 {E* , t*) given 
by Eq. (11) into Eq. (AlO) with E* = Tqx, finally we obtain the kinematical 
relation between 7 and n° pseudo-rapidity distributions 



drj* 



l + V^V^tanhN* 
2 \ Po ) F q (t v .) 



dN 7 
drj* 



(All) 



with 



A( V >) = / - I In f^21 - B„, v . ) (1 - x)>e-V 




■-;}' (A12) 

see Eq. (19) for r*. In Fig. (19), we present the numerical value of the 
correction rate for the Sternheimer approximation with a* w 0, corresponding 
to the second term in the square bracket in Eq. (All), where we assume po = 
200MeV/c, corresponding to approximately p^ = 185, 189, 191, 193MeV/c for 
v / s = 0.1, 0.2, 0.5, l.OTeV respectivey. We find that it is as large as 7% around 
j]* w 2 at Vs = 1 TeV. 
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